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' A perturbative treatment of reduced density operators of quantum subsystems is implemented 

, in the same spirit as Fermi Golden Rule for scattering. Analytic expressions for linear entropy (a 

. measure of purity loss, and in some cases of coherence loss) and for subsystem's energy variations 

04 ' (dissipation) are given. Application to electromagnetic field superposition states in a dissipative 

^ ^ cavity is performed. Evaluation of typical field and reservoir time scales (tf and tu) show that 

. even for small temperatures they are very different. We also indicate the condition on the cutoff 

temperature above which the decoupling assumption is quantitatively justified. 
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' Since 1970 many techniques for controlling single quantum systems have been developed, e.g., trapping a single 
^ ion in what is called an "ion trap" [ij, or storing a single mode of electromagnetic field in high quality cavities 13], 
trying to isolate them from the rest of the world and thus allowing us to probe many different aspects of its behavior 
. with amazing precision. However, the control over properties of quantum systems, the essential ingredient for the 
^\ ' technological implementation of quantum computers (as e.g. entanglement among many qubits), is strongly hindered 
. by deleterious environmental effects. The first experiment which clearly exhibits the effect of the environmental 
nI ' dynamics on a mesoscopic quantum superposition has been performed at the Laboratoire Kastler-Brosscl in Paris . 

From the theoretical point of view several independent ideas have been developed to treat open quantum systems 
The best known and used is the approach with master equations. Its starting point consists of a system of interest 
P-i. {S) Hnearly coupled to a "reservoir" (i?) usually simulated by a large number of harmonic oscillators (for a complete 
derivation see la, IgI, 'TiI) . The trace over the reservoir degrees of freedom leads to an effective equation for the system of 
^ interest containing nonunitary operators, responsible for coherence and energy loss, where the reservoir appears only 
^ through some mean values, just like in Thermodynamics. Several approximations are necessary in order to obtain 
such linear, markovian equation for S. One approximation which needs physical support to rest on -Sj is the one which 
^ ' considers the reservoir density uncoupled from S for all times J^. It can not be rigorously true, since decoherence of 
S system can only come from correlations with R. As discussed by Cohen- Tannoudji and co-workers [^, this is easy 
to accept it provided the temperature is large enough so that the typical time scales are very different. This tends 
to make any previous correlation between R and S irrelevant after a short time. However, we know from Schmidt's 
decomposition that if we have any pure bipartite system, the evolution of correlations is exactly the same for both 
subsystems. In this case, a "cold" reservoir consisting of all oscillators in the ground state belongs to this category 
and therefore the argument breaks down. In spite of this fact, master equations usually work extremely well also for 
zero temperature 0. Why? This is one of the long standing questions in the area of open quantum systems. 

The objective of the present contribution is to work out the typical time scales for variation of coherence and energy 
in any two coupled quantum systems. When one of the systems is a reservoir we obtain a limiting temperature for 
the validity of the factorization assumption holds. For this purpose we next implement a perturbative treatment of 
reduced matrices of quantum systems in the same spirit as the Golden Rule for scattering. 

Consider a Dyson like expansion for reduced operators which obey the following Hamiltonian dynamics 

H = Hi+H2 + Hc^Ho + Hc (1) 

where Hq contains the autonomous evolution of the systems in question and Hq represents the coupling between 
them. We will obtain perturbative equations for the reduced densities 



P^ = Trj73(i) , 



(2) 
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where p{t) is the total density, i = 1,2, j = 1,2 and i 7^ j are subsystem's indexes. We get 
Mt) = p,{0) - i f dt'Tr, [Hc{t'),p{Q)] 



dt' f d<"Tr,[i/c(<'), [Hc(t"),p(0)]] +•■• 
Jo 

— (0) , (1) , (2) 

= p\> +p\' +p\' +■■■ , 

where Hc{t) = o^^ot ^^^-iHgt j^^j-^g interaction representation is appUed). Note the symmetry in the indexes 
which will allow us to study the short time dynamics of any of the two subsystems. Instead of proceeding formally 
and rewriting the above equation in terms of matrix elements, we prefer to illustrate the procedure for obtaining 
dissipation and decoherence time scales by mean of an example. Let = 1) 

Ho = Loa) a + ^njb]bj , (4) 

j 

which represents, e.g., the free Hamiltonian of a given cavity mode plus that of the environment simulated by oscillators 
as usual, oj, Vlj, a(at) and are the natural frequencies, annihilation (creation) operators of the mode and the 

reservoir, respectively. We assume that the coupling between the two systems is given by 

He ^Y.^Aa%, + ah]) , (5) 

j 

where 7^ are the coupling constants assumed to be real. 

Subsystems dissipation timescales: Let us define a subsystem energy by 

E,{t) ^Tr [H,p,{t)] (6) 

and assume that we initially have the state of the composed system given by 

where f3 — l/{kBT). Inserting Eq. © in we get, up to second order for the cavity mode labeled 1, 

E,it) ^ e["^ + E['\t) + E^^\t) + ■ ■ ■ , (8) 

where the upper indexes indicate the order of the approximation. We get 

= Tri[tjatapi(0)] = uj{a''a)o = Ei{0), (9) 
pt 

E[^\t) = -iuj / dt'TT{a^a[Hcit'),p{0)]} = 0, (10) 
Jo 

Ei^\t) = ^LO f dt' f dt"TT{a^a[Hc{t'),[H,{t"),pm]] 
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■^-^ I 7j sin Ajf ^ ^ 



(ata)o) , (11) 



where Aj = (fij — cj) /2, fij = (e^^^ — l) ^ is the average occupation number of the j-th reservoir state at temperature 
T, and (0^0)0 is the average initial excitation number of the subsystem 1. Note in the above equation the basic 
ingredients of Fermi's Golden Rule: each reservoir mode contributes to the system energy initially in a way proportional 
to its coupling constant and to the difference in occupation numbers between the mode and the system. As time goes 
by, the (sinAjt/Aj)^ function selects only those modes with frequency very close to lo. In scattering language, this 
is the "phase space" available for the process. 

In order to gain deeper physical insight into the energy and coherence loss processes, let us consider the subsystem 
1 initially in the so called even coherent state jlOl | with amplitude a, pi (0) = \ea) {ea\, where 

lea) = r-p^ Wo') + I ^ ")] ■ 

' ' 2(1 + e-l"l') LI ' I 
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Here, |a) is Glauber's coherent state fllj. The average initial excitation number is given by 

(a^a)o = |ap tanh |ap . 

To proceed with the calculation we approximate (as usual) the summation over j by an integral in Eq. 



(12) 



sin Ajt 



dflg{n)F{fl) 

-l/(2t) 



sin At 



uj-l/(2t) 



sin At 



l/(2t) ^ 

d A g{uj + /\)F{uj + A) = — 

-l/(2t) '''F 



(13) 



where g(f2) is the density of modes of the reservoir at frequency and A = {U, ^ uj) /2. The first approximation 
is like a stationary phase argument and the second comes from a linear approximation to the sine function close to 
A = 0. This equation defines the time scale Tp^ provided F has frequency squared dimension. 
For the expansion in energy of subsystem 1, we get 



[t) ~ E, (0) 



cut 



7(f7) 



ujt''\a\"t'dnh\a\" / dng {n) 



sin {At) 



At 



sin (At) 
At 



(14) 



Note that in the above equation, the second term on the r.h.s. is independent of temperature and therefore defines a 
dissipation time scale, whereas the first term, on contrary, depends on temperature and allows us to define a thermal 
timescale. These timescales are determined through the ratio \Ei (t) — Ei (0)| /Ei (0), and reads 



t 



l/(2t) 



dA g {uj + A) 7^ {uj + A) 



(15) 



-1/(24) 



'th 



\af tanh \af 7-i/(2t) 



l/(2t) 



dA g (to + A) n{LU + A) 7^ {lu + A) . 



In the asymptotic case t — > cx), which corresponds to the Markovian approximation, we obtain 



and 



which give us the general timescale relation 



'th 



Tth 
Tdis 



= g(w)7^(w) 

n{uj)g{uj)j'^{u}) 
lal^ tanh lal^ 



lap tanh laP 



n (lu) 



(16) 



(17) 



(18) 



(19) 



In the large temperature regime f3u! ^1 we have Tth /"Tdis PoJ tanh |a|^ and Eq. ifT^I) has a simple semiclassical 
interpretation: as |a|^tanh|a|^ is proportional to oscillator's mean energy, the larger \a\^ is, the faster dissipation 
occurs; accordingly, as n (uj) = exp (— /3w) / [1 — exp (— /^w)] gives the mean number of quanta of frequency oj on 
thermal reservoir at temperature T, the larger n(a;), the faster thermal effects appear. Zero temperature is thus 
characterized by the absence of thermal effects. 

Subsystems decoherence timescales: The expression for the idempotency deficit |li] of the subsystem up to second 
order, is given by 



(t) = 1 - Tr.pi it) « it) + (t) + Sr' it) 



c(2) 



where ,5^^ (t) = ,5,(0), Sl'> (<) 



(1) 



-2Tr, 



last expressions, the operator p^"-* (t) represents the n-th term of the expansion For the state given by Eq. 0, 



p,(0)pf^ {t)],andSP (t) 



Tr, 



pf^ it) 



2Tr, 



P. (o)pr it) 



(20) 



In the 
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we have 5^^^ (t) — and J^"* (t) = — 2Tri pi (0) p^^^^* (t) . Assuming that the subsystem 1 was prepared in the initial 
state pi (0) — \ea) (ea|, the evaluation of S\ {t) yields, in the continuum limit, 



Si (t) 



dng [fl) 



sin (At) 



|n(f}) (^|a|^tanh|Q;|^ + 1 



[n(f^) + l]|a|^tanh|a|^}. (21) 



We can also employ the same reasoning that allows us to determinate the dissipation timescales in order to define the 
timescale t^cc- In the present case, purity loss and decoherence timescales coincide. Therefore, within the Markovian 
approximation, we get 



ll = 2g (iu) 7^ {uj) |n (uj) (^\a\^ tanh \a\^ + l) + [n (lj) + 1] \a\^ tanh 



(22) 



As is clear, in the expressions (|14() and (|21() the ingredients are the same, but the results are quite distinct. The 
greater difference is the relative sign and it is easy to understand it in this example. In both cases, the first term is the 
temperature contribution and the second is still there even for zero temperature. The first term tends to increase Ei 
(it is a "hot" term), while the second to decrease it (it is a "cold" one), but both add to decoherence. This difference 
in sign reflects in the difference in behavior and timescales of these processes. 

In case T = only "cold" terms contribute, and we obtain the relation 



2 lap tanh lap 



(23) 



and there is no thermalization in the sense which we are using this word. In the limit |ap ^ 1, Eq. (|23(l clearly shows 
the dependence of the purity loss timescale with the separation of the states that form the superposition in the initial 
state pi (0). In fact, the more distinguishable these states, the smaller the purity loss timescale. Such relations of 
time scales are usually considered the answer to the impossibility of naturally finding the great majority of quantum 
states for macroscopic systems (e.g.: Schrodinger cat states) [l3|. 
For temperature T we have 



Tth 

^dec 



2 lap tanh |ap 



n (u;) 



I [n {lo) + 1] |ap tanh |ap + n (uj) |^|ap tanh |ap + 



which stresses the fact that decoherence is even faster than thermalization. In the same way 



^ = 2 |[n (c^) + 1] lap tanh |ap + n (lu) (|ap tanh |ap + l) } , 



(24) 



(25) 



where |a| plays its crucial role, and for large temperatures we obtain the expected behavior of linear increasing of 
decoherence rate with respect to temperature, in complete agreement with expression (20) of ^14] and also with |l5| . 
but obtained in a very different framework. To authors knowledge, it is the first time general timescale relations as 
Eqs. (|24|l and H25|) appears for this important example. 

Master equation ~ a quantitative indicative of separability: The time scale of the reservoir can be straightforwardly 
evaluated by using the expansion Eq. (|2()|l now for the linear entropy of the reservoir S2{t) = 1 — Tr2/C2(i)^- If 
have Tr2[/92(0) — P2(^)]/Tr2/O2(0) ~ 0, we can consider that p2{t) « P2(0). However, up to second order we get 



^^2pl{t) 

Tr2pi(0) 



N 



^2 sin^ ( Afet) 



k=l 



{Akty 



4(2(ata) + l) 



tanh 



- 1 



+ 2 (a^a) - (a) (a^) tanh 



Using the same procedure as before one can obtain the characteristic coherence change for the reservoir as 



5(^)7^(^) 



'2. do 



2n 



/[^ ' |2 (atfl) [2n{uj) + 1] - 8n(w) {2{a^a) + I) - {a){a^)\ 
[uj) + 1 



A reasonable condition for the separability of the timescales of systems 1 and 2 is 

T2,dcc > , 



(26) 



(27) 



(28) 
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where u ^ is the characteristic time of the free evolution of system 1 . The worst case we have is when the temperature 
is zero and therefore n{Lo) = 0. In this regime we can guarantee the use of the approximation p2{t) ~ P2(0) if we have 

g{u;)-f^{ij) {2{a^a} - {a){a^)) < uj. (29) 

Just to have some magnitudes to work with, we took from Ref. (sj the values {a} a) sa 9.5, n(ijj) — .05, tojl-K = 51 GHz, 
and Tdis = 160 /^s, and obtain T2,dcc « 13 /iS ^ u)^^ « 20 ps, as one would expect. 

If we want to have a condition which allows for the determination of a cutoff temperature in terms of dynamical 
ingredients and initial conditions we should use the condition (|28|l in Eq. (|27|l . 

To summarize we have provided a quantitative support for the factorization assumption between the system of 
interest and the reservoir density matrices which is widely used in the decoherence literature. Moreover we give 
analytical expressions for typical correlation and relaxation timescales which depend only on general ingredients such 
as the Hamiltonian and the initial state. 
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